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The specific behavior of spinor field in curve space-time with the exception of FRW model 
almost always gives rise to non-trivial non-diagonal components of the energy-momentum 
tensor. This non-triviality of non-diagonal components of the energy-momentum tensor 
imposes some severe restrictions either on the spinor field or on the metric functions. In this 
paper within the scope of an anisotropic Bianchi type-I Universe we study the role of spinor 
field in the evolution of the Universe. It is found that there exist two possibilities. In one 
. scenario the initially anisotropic Universe evolves into an isotropic one asymptotically, but 

["t | ! in this case the spinor field itself undergoes some severe restrictions. In the second scenario 

\ the isotropization takes places almost at the beginning of the process. 
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I. INTRODUCTION 

According to the inflationary scenario, it is believed that a scalar field known as inflaton is 
responsible for a rapid accelerated expansion of the early Universe [1-31. For the inflationary 
mechanism to work there must exist a weakly coupled scalar field Wicp is initially at a false 
vacuum which leads to the inflation until the right vacuum value is (^t|a|ned. This inflationary 
model solves the problem of flatness, isotropy of microwave background radiation and unwanted 
relics. Contrary to the prediction of the standard cosmological models, recent observations showed 
an accelerated mode of expansion of the present day Universe [4, 5]. Though the existence of an 
inflationary scenario is not of much concern, the question of whers the scalar field comes from 
and why it undergoes such a peculiar phase transition from fa!}^ right vacuum still remains 
unanswered. This leads cosmologists to reconsider alternative possibilities. 

As one of the way out many specialists considered spinor field as an alternative source. Being 
related to almost all stable elementary particles such as proton, electron and neutrino, spinor field, 
especially Dirac spin-1 /2 play a principal role at the microlevel. However, in cosmology, the role 
of spinor field was generally considered to be restricted. Only recently, after some remarkable 
works by different authors [6-20], showing the important role that spinor fields play on the evo- 
lution of the Universe, the situation began to change. This change of attitude is directly related to 
some fundamental question|_pyriodern cosmology: (i) problem of initial singularity; (ii) problem 
of isotropization and (iii) late time acceleration of the Universe. 

(i) Problem of initial singularity: One of the problems of modern cosmology is the presence 
of initial singularity, which means the finiteness of time. The main purpose of introducing a 
nonlinear term in the spinor field Lagrangian is to study the possibility of the elimination of initial 
singularity. In a number of papers r8— 121 it was shown that the introduction of spinor field with a 
suitable nonlinearity into the system indeed gives rise to singularity-free models of the Universe. 

(ii) problem of isotropization: Although the Universe seems homogenous and isotropic at 
present, it does not necessarily mean that it is also suitable for a description of the early stages of 
the development of the Universe and there are no observational data guaranteeing the isotropy in 
the era prior to the recombination. In fact, there are theoretical arguments that support the existence 
of an anisotropic phase that approaches an isotropic one [211 . The observations from Cosmic 
Background Explorer's differential radiometer have detected and measured cosmic microwave 
background anisotropics in different angular scales. Thes$ apisotropies are supposed to hide in 
their fold the entire history of cosmic evolution dating back to the recombination era and are being 
considered as indicative of the geometry and the content of the universe. More about cosmic 
microwave background anisotropy is expected to be uncovered by the investigations of microwave 
anisotropy probe. There is widespread consensus among the cosmologists that cosmic microwave 
background anisotropics in small angular scales have the key to the formation of discrete structure. 
It was found that the introduction of nonlinear spinor field accelerates the isotropization process 
of the initially anisotropic Universe [ 10. 11. 131 . 
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(iii) late time acceleration of the Un iverse: Some recent experiments detected an accelerated 
mode of expansion of the Universe 1 4 , 111. 



;ection and further experimental reconfirmation of 
current cosmic acceleration pose to cMrhblogy a fundamental task of identifying and revealing the 
cause of such phenomenon. This fa<gt| be reconciled with the theory if one assumes that the 
Universe id mostly filled with so-called dark energy. This form of matter (energy) is not observable 
in laboratory and it does not interact with electromagnetic radiation. These facts played decisive 
role in naming this object. In contrast to dark matter, dark energy is uniformly distributed over the 
space, does not intertwine under the influence of gravity in all scales and it has a strong negative 
pressure of the order of energy density. Based on these properties, cosmologists have suggested a 
number of dark energy models those are able to explain the current accelerated phase of expansion 
of the Universe. In this connection a series of papers appeared recently in the literature, where a 
spinor field was considered as an alternative model for dark energy [14 —16. 181 . 

It should be noted that most of the works mentioned above were qarffied c ut within the scope of 
Bianchi type-I cosmological model. Results obtained using a spinoy fjqld isji source of Bianchi 
type-I cosmological field can be summed up as follows: A suitable choice of spinor field nonlin- 
earity 

(i) accelerates the isotropization process [10. 11. 131: 
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(ii) gives rise to a singularity-free Universe [10-131; 

(iii) generates late time acceleration [14-16] laj 191. 

Given the role that spinor field can play in thm evolution of the Universe, question that naturally 
pops up is, if the spinor field can redraw |t ihe|pic| t ure| of evolution caused by perfect fluid and 
dark energy, is it possible to simulate perfect fluid and dark energy by means of a spinor field? 
Affirmative answer to this question was given in the a number of papers ["22-261 . In those papers 
spinor description of matter such as perfect fluid and dark energy was givenaM tHe evolution of the 
Universe given by different Bianchi models was thoroughly studied. In almost all the papers the 
spinor field was considered to be time-dependent functions and its energy-momentum tensor was 
given by the diagonal elements only. Some latest study shows that due to the specific connection 
with gravitational field the energy-momentum tensor of the spinor field possesses non-trivial non- 
diagonal components as well. In this paper we study the role of non-diagonal components of the 
energy-momentum tensor of the spinor field in the evolution of the Universe. To our knowledge 
such study was never done previously. In section II we give the spinor field Lagrangian in details. 
In section III the system of Einstein-Dirac equations is solved for BI metric without engaging the 
non-diagonal components of energy-momentum tensor as it was done in previous works of many 
authors. In section IV we analyze the role of non-diagonal components of energy-momentum 
tensor on the evolution of the Universe. 



H. SPINOR FIELD LAGRANGIAN 

For a spinor field l//, the symmetry between \y and iff appears to demand that one should choose 
the symmetrized Lagrangian [27]. Keeping this in mind we choose the spinor field Lagrangian as 



[10] 



m sp v/v/-F, (2.1) 



where the nonlinear term F describes the self-interaction of a spinor field and can be presented 
as some arbitrary functions of invariants generated from the real bilinear forms of a spinor field. 
Since i// and yr* (complex conjugate of \jf) have four component each, one can construct 4x4= 16 
independent bilinear combinations. They are 

S = \[r\if (scalar), (2.2a) 

P = i\j/y 5 \j/ (pseudoscalar), (2.2b) 

= (v/y'V) (vector), (2.2c) 

= (y^y^y) (pseudovector) , (2.2d) 

T^ v = (xj/o^ V) (antisymmetric tensor), (2.2e) 

where o^ v = (j/2)[y^y v — y v y"]. Invariants, corresponding to the bilinear forms, are 

(2.3a) 
(2.3b) 
(2.3c) 

(2.3d) 

(2.3e) 

According to the Fierz identity, among the five invariants only / and J are independent as 
all others can be expressed by them: I v = —Ia = I + J and 7j = I — J. Therefore, we choose 
the nonlinear term F to be the function of / and J only, i.e., F = F{I,J), thus claiming that it 



/ 


= s\ 




J 


= P\ 




Iv 


= v^ = 


(Wfv) 8tiv(Y7 v ¥), 


Ia 


= A M A" = 




h 


= W v 


= (yo^^g^agvpiYO^V) 



4 



Saha B. 



describes the nonlinearity in its most general form. Indeed, without losing generality we can 
choose F =F(K), with K={I,J,I + J,I — J}. Here is the covariant derivative of spinor field: 



dxif d\lf 

v MV = ^r-r M v, v^^+^rv, (2.4) 



with r^x being the spinor affine connection. In (12.11 ) y's are the Dirac matrices in curve space-time 
and obey the following algebra 

f r v + r vf =2g ^ (2.5) 
and are connected with the flat space-time Dirac matrices y in the following way 

8pv(x) = el(x)e b v (x)r\ ab , y M (*) = ej(x)y fl , ( 2 -6) 

where r\ ab = diag(l, — 1,-1,-1) and is a set of tetrad 4- vectors. The spinor affine connec- 
tion matrices (x) are uniquely determined up to an additive multiple of the unit matrix by the 
equation 

V M y v = -T%y p -T M y v + y v T fl = 0, (2.7) 

with the solution 

r^lfafd^-^ypfT^ (2.8) 

Varying (12.11 ) with respect to vKvO one finds the spinor field equations: 

iy^V^xj/- m sp xj/ - 2F K (SK! + iPKjy 5 )^ = 0, (2.9a) 
+m sp x[r + 2F K x[f(SK 1 + iPKjy 5 ) =0. (2.9b) 

Here we denote Fk = dF /dK, Kj = dK/dl and Kj = dK/dJ. 

We consider the case when the spinor field depends on t only. In this case for the components 
of energy-momentum tensor we find 

T^ = m, v S+F(K), (2.10a) 
T} = T 2 2 = T 3 3 =F(K)- 2KF K . (2. 10b) 

The energy-momentum tensor of the spinor field is given by 

T$ = (nVvY+WvV»Y-V»W V Y-VvWiiA -SjLsp (2.11) 
where L sp in view of (12.91 ) can be rewritten as 

Lsp = ^[¥Y^^¥-^^¥f¥] -m sp \frY-F(K) 

= ~ Wlt^ii W-ms V v}- 1 - [V^ xjrf + m sp \j/} yf—F (K) , 

= 2F K {IKi + JKj) —F = 2KF K —F(K). (2.12) 

Let us now recall that in the unified nonlinear spinor theory of Heisenberg, the massive term 
remains absent, and according to Heisenberg, the particle mass should be obtained as a result 
of quantization of spinor prematter [28. 291 . In the nonlinear generalization of classical field 
equations, the massive term does not possess the significance that it possesses in the linear one, 
as it by no means defines total energyf (o|r|mftss) of the nonlinea r fiel d system. Moreover, it was 
established that only a massless spinor field with the Lagrangian (12.11) describes perfect fluid from 
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phantom to ekpyrotic matter [22-261. Thus without losing the generality we can consider the 
massless spi nor fiel d putt ing tn = ). 

Inserting (12. 10a!) and (12.10q|) ip|to j he barotropic equation of state 



where W is a constant, one finds 



p = We, 



2KF K =(1+W)F(K), 



(2.13) 



(2.14) 



with the solution 

F(K) = A^ 1+H/ )/ 2 , A = const. (2.15) 

Depending on the value of W (12.131) describes perfect fluid from phantom to ekpyrotic matter, 
namely 



W = Q, 


(dust), 


(2.16a) 


W = l/3, 


(radiation) , 


(2.16b) 


We (1/3,1), 


(hard Universe), 


(2.16c) 


W = l, 


(stiff matter), 


(2.16d) 


We (-1/3,-1), 


(quintessence), 


(2.16e) 


W = -l, 


(cosmological constant), 


(2.16f) 


W<-1, 


(phantom matter) , 


(2.16g) 


W>1, 


(ekpyrotic matter). 


(2.16h) 



In account of it the spinor field Lagrangian now reads 



A ^(lW)/2_ 



(2.17) 



Thus a massless spinor field with the Lagrangian (12.171 ) describes perfect fluid from phantom to 
ekpyrotic matter. Here the constant of integration A can be viewed as constant of self-coupling. A 
detailed analysis of this study was given in [22-251. 

A Chaplygin gas is usually described by a ec uation of state 



p = -A/ 

Then in case of a massless spinor field for F one finds 

FUF __ 1 dK 
F l +r-A ~ 2~K' 



with the solution [23 



251 



F={A + XK^ 2 ) mi+r) . 



The Spinor field Lagrangian in this case takes the form 



(A + A^ 1+ ^/ 2 ) 1/(1+r) . 



^~ 2 

Finally, it should be noted that a quintessence with a modified equation of state 



(2.18) 



(2.19) 



(2.20) 



(2.21) 



p = w(e-£ CI ), we (-1,0), 



(2.22) 



6 



Saha B. 



where £ cr some critical energy density, the spinor field nonlinearity takes the form 

W 



F = Xs i+w + 



The spinor field Lagrangian in this case reads 



l+W 



L 



sp 



XK (i+w)/2_. 



W 



(2.23) 



(2.24) 



Setting e cr = one gets (12.171) . The purpose of introducing the modified EoS was to avoid the 
problem of eternal acceleration. 

A detailed study of nonlinear spinor field was carried out in [10 . 11, 13J. In what follows, 
exploiting the equation of states we find the concrete form of F wliic i 
perfect fluid and dark energy. 



des 



ri ?es various types of 



III. BIANCHI TYPE-I ANISOTROPIC COSMOLOGICAL MODEL 

Let us study the evolution of the Universe filled with spinor field. In doing so we consider the 
case when the gravitational field is given by an anisotropic Bianchi type-I cosmological model. 

Bianchi type-I (BI) model is the simplest anisotropic cosmological model and gives an excellent 
scope to take into account the initial anisotropy of the Universe. Given the importance of BI model 
to study the effects of initial anisotropy in the evolution of he Universe, we study this models in 
details. 

We consider the BI metric in the form 



ds 2 = dt 2 — a 2 dx 2 — a\ dy 2 — dz 2 , 

with a, b and c being the functions of time only. 

For further purpose we define the volume scale V of the BI metric as 

V = abc. 

The system of Einstein equations in this case reads 



(3.1) 



(3.2) 



a 2 a 3 ^_ a 2 a 3 


= KTi, 


(3.3a) 


Cl2 CI3 &2 ®3 


#3 d\ ^dj, d\ 


= kT 2 , 


(3.3b) 


a?, a\ aj, a\ 






(X\ &2 CL\ d2 


= KZ3 3 , 


(3.3c) 


a 1 (12 Ci\ d2 






®1 @2 ^_@2@3 @3 ®1 


= <• 


(3.3d) 


Q\ d2 Cl2 0,3 CI3 (l\ 





Here is the energy momentum tensor of the spinor field. 

Solving the Einstein equation on account of the fact that = T 2 = for the metric functions 
one finds [101 



Cl; 



AV 1/3 exp(W^), nA = l, E* = 0, 
v J v J i=l (=1 



(3.4) 



with Di and X ; being the integration constants. Thus we se that the metric functions can be ex- 
pressed in terms of V. 
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Summation of (I3.3al) . (I3.3bl) . (I3.3cl) and 3 times (I3.3dl) leads to the equation for V [101 



v = ^(r° +t})v. 



(3.5) 



As we have already found, the c omp onents of energy momentum tensor are the function of K. If K 
is a function of V, then the Eq. (|3.5I) possesses exact solu tion. In order to show that K is a function 
of V we go back to the spinor field equations. From (12.91 ) one dully finds 

S -4F K PKjA = 0, (3.6a) 
P + 4F K SK I A = 0, (3.6b) 
A - AFkSK^q + AF k PKjSq = 0, (3.6c) 

where Sq = SV, P = PV, A = AV with A = v/y°yV- Summation of (I3.6al) . (I3.6bl) and (I3.6cl) 
leads to 

S 2 +P 2 +A 2 =Ci/V 2 , Ci=const. (3.7) 
On the other hand from (|3.6al ) and (|3.6bl ) one finds 

£/S 5 + Wo = 0. (3.8) 
In case of K = I, i.e., Kj = 1 and Kj = from (13.81 ) one finds 

K = I = S 2 =C//V 2 , C/ = consf. (3.9) 
For K = J, i.e., = and = 1 from (13.81 ) one finds 

K = J = P 2 = Cj/V 2 , Cj = const. (3.10) 
If if = / + 7, i.e., = 1 and Kj = 1 from (HU) one finds 

K = I + J = S 2 + P 2 = C I+J /V 2 1 C 1+J = const. (3.11) 
and finally, for K = I — J, i.e., Kj = 1 and A"/ = — 1 from (13.81) one finds 

K = I-J = S 2 -P 2 = Q j/V 2 , Q-j = const. (3.12) 
Thus we see that for the BI spacetime given by (13.11 ) one finds 

V 2 

K=A, V = const. (3.13) 
V z 

In case of (|2.17l) we have 

r ° = e = A^ 1+M/ )/ 2 , (3.14a) 

T { 1 = -p = ~We = -WXK^ l+w ^ 2 . (3.14b) 



Eq. (13.51) then takes the form 
with the solution in quadrature 



V = ^XV ( l +w (l-W)V- w , (3.15) 



dX t+tQ. (3.16) 



3kW* +w V 1 - w +d 
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Here C\ and to are the integration constants. 

Let us consider the case when the spinor field is given by the Lagrangian (12.211) . In this case 
we have 



T ° = e=(A + XK^y 2 ) W+7 \ 
T\ =- p =A/ £ v=A/(A + XK^+y)' 2 ) r/(1+r) . 
The equation for V now reads 



V 



3k 



(AV l+ v+XV l+r ) mi+r) +AV l+ y(AV l+ v + XV l+Y ) 



i+y\y/(i+y) 



(3.17a) 
(3.17b) 



(3.18) 



with the solution 



dV 



Ci +3KV(AV l +r + XVq 7 ) 



i+ni/(i+y) 



t + to, C] = const, to = const. 



(3.19) 



Inserting y = 1 we come to the result obtained in DPI . 

Finally we consider the case with modified quintessence. Taking into account that 

r = U (W)/2 + JL £cr) 
o i +w cr ' 



r/ = r 2 2 = r 3 3 



+ 



w 



l+W 



■'CIl 



for V in this case we find 
with the solution in quadrature 



XV Q 1+W (\-W)V- W +2We CI V/(l+W) 



-w 



I 



dV 



3K[lvi- w V l - w +We a V 2 /(l+W)]+Ci 



t+tQ. 



(3.20a) 
(3.20b) 

(3.21) 
(3.22) 



Here C\ and to are the integration constants. Comparing (|3.22l ) with those with a negative A-term 
we see that £ cr plays the role of a negative cosmological constant. 

Let us also write the components of the spinor field explicitly. For simplicity, we consider the 
case when K = I. Defining <2i = 2SFkKj = 2SFk in this case we find 



Xj/^t) = (Ci/Vy)exp y-ij 9dt 
Wiit) = (C 2 /W)exp i^-ij @dt\, 
V3(0 = (C 3 /VV)exp (ij ®dt 
y^(t) = (C 4 /W)exp I ij 9dt 



(3.23a) 
(3.23b) 
(3.23c) 
(3.23d) 
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with Ci, C2, C3, C4 being the integration constants and related to Vq as 

! Ci +C 2 C2 — L3C3 — C4C4 = Vo- 

In what follows we will study the obtained results within the scope of some recent findings, 
namely the fact that the spinor field possesses noz-trivial non-diagonal components of the energy- 
momentum tensor. 



IV. WHAT'S WRONG? 

In first view everything looks good and the papers written till the date on this subject seems 
correct. But there is still something to be worried about. In what follows, we speak about the new 
findings on this field. 

It should be remembered that the spinor field is more sensitive to the gravitational one. It is due 
to specific spinor connection in curve space-time. So, let us first write the spin affine connection 
explicitly. For BI metric it looks: 

r = o, r i = ^ff 1 (4.i) 

a 

Taking it into account from (12.1 II) it can be easily verified that the energy-momentum tensor of the 
spinor field possesses non-trivial non-diagonal components as well. 

T^ = m sp S + F(K)=F(K), (4.2a) 
T} = T 2 2 = T 3 3 = 2KF K —F(K), (4.2b) 

r 2 ' = p(^-%rWV, (4.20 

T] = ^-lh- a J.) f fff V , (4.2d) 
T} = ^h-^)yfffy. (4.2e) 

4 «2 \Cl2 a?, J 

So the complete set of Einstein equation for BI metric should be 

* + * + *l*l = K{ ZKF K -F{K)), (4.3a) 

a,2 C13 ci2 0,3 

=fc(2HW(Z)), (4.3b) 

CI3 a.\ $3 a,\ 

^ + ^ + ^ = K (2KF K -F(K)), (4.3c) 

Ci\ a.2 Cl\<X2 



(l\ 0,2 Cl2 a 3 a 3 a \ 

I CL2 ( a\ 0.2 



0=pp-- WWV, (4-3e) 
4 a\ \a\ a.2 J 

4 a\ \a3 a\ J 

4 «2 \ a 2 ®3 / 
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The equations (I4.3el) . (I4.3fl) and (|4.3g ) impose some severe restrictions either on the spinor field, 
or on the metric functions, or on both of them. 

If the restrictions are imposed on the spinor field, we obtain 



yy 1 y 2 /V = v / 7 3 y 1 y ( V = v A 7 2 7 3 y°v / = o. 



(4.4) 



In this case the expressions obtained for the metric functions in the earlier papers remain unaf- 
fected. But the components of the spinor fi eld w ill undergo some changes. It should be verified 
that in this case the integration constants in (|3.23l ) should obey 



lM — C 2 C2 — C3C3 +C4C4 

C\ C 2 + Cld + C 3 *C 4 + C4C3 

C*/— 1 1 /'"<*/'"< S"'*/^ s^*/^ 

|Ci +02^2 — C3C3 — C4C4 



0, 
0, 
0, 

Vb, 



which gives 



C1C2 + C3C4 



CIQ+C4C3 



2C2 — C4C4 = — 



0, 
2 



(4.5a) 
(4.5b) 
(4.5c) 
(4.5d) 



(4.6a) 
(4.6b) 



I n the c ases considered above for volume scale we obtained the expressions given by (I3.16I) , (I3.19I ) 
and (|3.22l) . for the Universe filled with quintessence, Chaplygin gas and quintessence with mod- 
ified equation of state, respectively. The Universe in these cases is initially anisotropic which 
evolves into an isotropic one asymptotically [ 25. 261 . 

}f the spinor field unaltered. In this case from 



The other poss ibility is to keep the components 



(I4.3el) . (I4.3fl) and ( |4.3g| for the metric functions 



qmq immediately finds: 



a\ ci2 a-i _ a 
a\ a,2 az a 



Taking into account that 



a. 



d fdi 



cii dt\aj/ w// 



+ 



(i;\ 2 d /a\ /d < 2 
dt\aJ \a< 



the system (14.31 ) can be written as a system of two equations: 



„a a 
2- + - 

a a z 



1 ! 



a 

* 

a 



(4.7) 



(4.8a) 
(4.8b) 



In order to find the solution that satisfies both (I4.8al ) and (I4.8bl) we rewrite (I4.8al > in view of 
(I4.8bl) in the following form: 

a = ^(3T} -T$)a. (4.9) 



Thus in account of non-diagonal components of the spinor field, we though begin with Bianchi 
type-I space time , in reality solving the Einstein field equations for FRW model. Before solving 
the equation (14.91) , let us go back to (13.41) . Taking into account that 
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in view of (14.71) we find that 

X l= X 2 =X 3 =0. (4.11) 

The trivi ality of the integration constant Xi follows from the fact that X\ + X 2 + X3 = 0. Thus the 
solution (13.41) should be written as 

3 

a i =D i V 1 / 3 =D i fi, Y[ D=1 ( 4 - 12 ) 

i=l 

which means it represents a tiny sector of the general solutions (|3.4I) which one obtains for the BI 
model in case of isotropic distribution of matter with trivial non-diagonal components of energy- 
momentum tensor, e.g., when the Universe is filled with perfect fluid, dark energy etc. 

Let us now define a = V 1 / 3 for different cases. In doing so we recall that K in this case takes 
the form 

K = -%, a = const.. (4.13) 
a b 

Then then equation for a in case of the spinor field given by (|2.17l) takes the form 

< 

with the solution is quadrature 



a = -^{\+W)af +W) a-( 2+m \ (4.14) 



{KX/3)af +W) a-^)+E x 

with E\ being integration constant. 

As far as Chaplygin scenario is concerned in this case we have 

k 2Aa 3 d+r) - Xaf +Y) 



da 

t, (4.15) 



6 a 2 (Afl3(i+r)-Aflf +r) ) 



3(1+7) \ 7/(1+7)" 



(4.16) 



This equation can be solved numerically. 

Finally we consider the case with modified quintessence. Inserting (I3.20al) and (I3.20bl ) into 
(14.91) in this case we find 



K 

d = ~6 



{ 3W + l)Xal^a-^-^e cr a 



with he solution 



+ 

da 



(4.17) 



(jc/3) \Xa 3 {l+W) a-(i w +V + [W/(l + W)]e ct a 2 + E 2 



t, E 2 = const. (4.18) 



It can be shown that in case of modified quintessence the pressure is sign alternating. As a result 
we have a cyclic mode of evolution. 

In what follows we illustrate the evolution of the Universe filled with quintessence, Chaplygin 
gas and quintessence with modified equation of state for two different cases: when the restric- 
tions are imposed on the spinor field and when the metric functions were restricted. In Figures 
[Q [2] and [3] we illustrated the evolution of the Universe filled with quintessence, Chaplygin gas 
and quintessence with modified equation of state, respectively. The solid (red) line stands for the 
volume scale, when the restrictions due to non-zero non-diagonal components of the energy mo- 
mentum tensor of the spinor field, were imposed on the components of the spinor field. In this case 
the isotropization takes place asymptotically. The blue line shows the evolution of the Universe 
when due to the non-zero non-diagonal components of the energy momentum tensor of the spinor 
field leads to the immediate isotropization of the Universe. Here we plot the volume scale as a 3 , 
which a being the average scale factor. 

As one sees, in case of early isotropization the Universe grows rapidly. 
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FIG. 1. Evolution of the Universe filled with quintessence. The solid (red) line stands for volume scale V, 
while the dash-dot (blue) line stands for a 3 . 




FIG. 2. Evolution of the Universe filled with Chaplygin gas. The solid (red) line stands for volume scale V, 
while the dash-dot (blue) line stands for a 3 . 



V. CONCLUSION 

Within the scope of Bianchi type-I space time we study the role of spinor field on the evolu- 
tion of the Universe. It is shown that even in case of space independent of the spinor field it still 
possesses non-zero non-diagonal components of energy-momentum tensor thanks to its specific 
relation with gravitational field. This fact plays vital role on the evolution of the Universe. There 
might be two different scenarios. In one case only the components of the spinor field are affected 
leaving the space-time initially anisotropic that evolves into an isotropic one asymptotically. Ac- 
cording to the second scenario, the space-time becomes isotropic right from the beginning,i.e., 



(5.1) 



Nonlinear Spinor Fields in Bianchi type-I spacetime reexamined 
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FIG. 3. Evolution of the Universe rilled with quintessence with modified equation of state. The solid (red) 
line stands for volume scale V, while the dash-dot (blue) line stands for a 3 . 



and can be completely described by the Einstein field equations for FRW metric. As numerical 
analysis shows, in case of early isotropization the Universe expands rather rapidly. There might 
be another possibility when the non-diagonal components of energy-momentum tensor influence 
both the spinor field and metric functions simultaneously. We plan to study this last case within 
the scope of other Bianchi models in our forthcoming papers. 
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